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^ ■ By using the example of photon-photon scattering box diagrams, different rep- 

resentations for loop integrals are discussed. A connection between hypergeometric 
^ , representation and dilogarithms is derived. An explicit analytic continuation for- 

^ ' mula to large values of Mandelstam variables is constructed. 

(N ' 

m 

^ ' 1. When studying many important physical problems (calculation of radiative corrections to 

' elementary particles interactions, examination of operator expansions, etc.), we are often con- 
fronted with various types of Feynman loop diagrams. The greatest success has been achieved 
in evaluating massless propagator-type diagrams (depending on one external momentum only). 
For such diagrams, this dependence is power-like, and the only problem is how to calculate 
numerical factors at these powers. Some powerful methods (Gegenbauer polynomial technique 
^ ■ m, integration by parts 0, the "uniqueness" method 0, and some others) have been developed 
to calculate these numbers, and for some classes of diagrams the results have been obtained for 
^ ! any numbers of loops (see, e.g., [^]). 

^ I The problem gets essentially more complicated, however, when we consider diagrams with 

masses and (or) with larger number of external lines. In this cases, dimensionless combinations 
can be constructed from squared external momenta and masses. As a result, we can get 
complicated functions of these variables (instead of numbers), and the calculation becomes 
non-trivial even on one-loop level. Thus, the problem is how to describe the occurring functions 
in all possible regions of momenta and masses. There are different approaches to this problem. 

Commonly-used "standard" representations for loop diagrams usually involve special func- 
tions like polylogarithms (LIat), generalized Clausen's functions, etc. Explicit results of such 
type are known for one- loop diagrams up to four external lines (see, e.g., in [§],§[) and for 
some particular cases of two- loop diagrams (see, e.g., 0, §). For some more complicated cases, 
however, only parametric integral representations are known (see, e.g., [0). 

There exists also another way to represent results for loop diagrams - in terms of hyper- 
geometric functions. In the convergence regions, these functions can be represented as simple 
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power series |TD|. Some results of such type have been presented in |n]. Moreover, in such a 
way it was possible to obtain explicit formulae for the general case of one-loop diagrams with 
arbitrary number of external lines |T^. Note that a special approach makes it also possible 
to represent integrals with tensor structure in numerators in the same form. If one need to 
know results in other regions of variables, it is possible to use analytic continuation formulae 
for occurring hypergeometric functions (see, e.g., in |14|). On the other hand, for complicated 
cases there is a problem how to construct analytic continuation to all regions of interest. 

To our opinion, both representations can be useful in realistic calculations. Because of this, 
it is also important to know how they are connected (how to go from one representation to 
another). Usually it is possible to use parametric integral representations for hypergeometric 
functions, and then to try to evaluate them in terms of known special functions (for example, 
polylogarithms). Such a procedure has been already considered by taking some examples of 
triangle diagrams [0, |l^ . In the present paper we would like to illustrate it via an example of 
box diagrams contributing to photon-photon scattering. 



2. Photon-photon scattering (77 — >■ 77) is the well-known example of non-linear effects 
in QED (see, e.g., - and references therein). Note that the processes gg — > 77 and 
99 ^ 99 s-re also described by the same diagrams. The lowest order of this process contains 
the box diagram presented in Fig. 1 as well as those obtained by permutations {k^ ^ k^) and 



Figure 1: 

(k2 <-> k^). Solid internal lines correspond to massive fermions (electrons or quarks) with the 
mass m, while dashed external lines denote massless particles (photons or gluons): 

—— ~~ ^3 ~~ ^4 ~~ 0* ("^) 

To simplify the analysis, in what follows we shall restrict ourselves by scalar integrals corre- 
sponding to these box diagrams, because it is known that integrals with tensor numerators can 
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be reduced to scalar ones (see, e.g., [2^, Formally, such scalar integrals correspond also to 
massless particles scattering in the Wick-Cutkosky model. 
As usual, we define the Mandelstam variables as 

s = {ki + k2f, t={k2-hf, u={ki-k^f- s + t + u = Q. (2) 

Note that only two of them (for example, s and t) are independent. So, if we denote the 
integral corresponding to Fig. 1 as J(s, t; m) then the whole amplitude (including diagrams 
with permutations) is proportional to the sum 

J{s,t]m) + J{s,U]m) + J{u,t;m). (3) 

Hypergeometric representation for J in four dimensions can be easily obtained from general 
formulae of the papers (see also |]TB|, ITDI) by putting kf = 0, 



,2 



t7T' 

J(s,t;m) = ^F3(l,l,l,l; 5/2 



t 



4m2 Am? 



(4) 



where 

F3(a..'.ft/3';.|..,).£E^ '"'-?'^'-'^''' (5) 



is Appell's hypergeometric function of two variables (see, e.g., ^^), and 



_ T{a + j) 

denotes the Pochhammer symbol. Note that the convergence region of the double series (|^) is 
X < 1, y < 1. So, formulae (§)-(||) yield asymptotic expansion of J{s,t;m) for small values of 
X = s/ {4m?) and y = t/ (4m^). 

To obtain "standard" representation for J{s,t] m), it is convenient to use the known para- 
metric integral representation for (see, e.g. 

r(7) 

Fs (a, a', /?,/?'; 7 1 x,y) - 



r(/3) Ti/3')Ti^- 13-/3') 



X / / d^dr] e^-^ 7]^'-^ (1 - e - v)^-''-^'-' (1 - ^x)-° (1 - w)""' (7) 



5>0, r;>0 



(it is valid for Re/5 > 0, Re/3' > 0, Re(7 - (3 - f]') > 0). For the case of §) we find, from (0), 
that 



F3 (1,1, 1,1; 5/2| x,t/) 



€>0, r;>0 



(i^ dr] 



(1 - ^a;)(l - r]y)^l-^-ri' 
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Evaluating the integrals on the r.h.s. of (^), we arrive at the following result: 
F3 (1,1, 1,1; h/2\x,y) 



+ E 

i=x,y 



2 In' 



2Li2 



+ In 



(^xy /3x 
Pxy ~l~ Px J 



In 



f^xy ~l~ /5y 



TT 

Y 



1 



-2Li2 



l^xy Pi 

A + 1 



In^ 



A + 1 



where 



Px^\l 



X 



y 



X y 



(9) 



(10) 



By simple transformations of dilogarithms, it is easy to show that the result (^), (^ coincides 
with the well-known expression (see, e.g., |]18|). Our expression, however, is more compact and 
contains four dilogarithms only (instead of eight ones in |]T8| ). 



3. To obtain the formula of analytic continuation to the variables 1/x and 1/y, it is conve- 
nient to use double Mellin-Barnes representation for F3, 



too too 



-too —too 



xr(-a) r(-6) r(« + a) T{a' + b) T{P + a) T{P' + b) ^ 

r(7 + a + b) 

where the integration contours are chosen so as to separate "right" and "left" series of poles 
of gamma functions in the integrand. If |a;| < 1, |y| < 1, we should close both contours to 
the right, and we obtain the formula (|^). On the other hand, for large and we should 
close the contours to the left. As a result, we obtain the following known analytic continuation 
formula (see, e.g., p^): 



{A,t^,a} r(a) r(7 - A - /i) v x 



XF2 yX + fx - 'J + 1, A, /i; A-p + 1, /i-(T + l 
where F2 is another Appell's hypergeometric function of two variables. 



1 1 

y. 



:i2) 



(a pp'.,,'luv)^ft^ (13) 



4 



and the sum (in ([T^)) extends over the following four sets of { A, /x, p, cr } : {a,a', P, P'}, 
{a, j3', j3, a'}, {P, a', a, /?'} and {/?, /?', a, a'}. Note that the radius of convergence of the function 
(^) isu + v <1. 

The case a = a' = P = P' = 1, 7 = 5/2 corresponds to so-called "logarithmic" case of 
analytic continuation, because in this case gamma functions in separate terms on the r.h.s. of 
([T^) have singularities (which cancel in the whole sum). Considering this limit of the expression 
([I^), we obtain the following analytic continuation formula: 



F3 (1,1, 1,1; 5/2\x,y) 



Axy 



ln(-4x) \n{-Ay) - — 



F2 1/2, 1,1; 1,1 



1 1 

x' y 



[ln(-4x) + \n{-Ay)] S^Fs - \n{-Ay) d^F2 - \n{-Ax) SyFs 



+dlF2 + do^d^F2 + do^d^>F2 + d^d^>F2 



(14) 



where we have introduced notation for parametric derivatives of the function ([13D , for example: 



d^F^^-^F^Up^P'; 7,7' 



1 


3 


X 





(15) 



a=l/2, /3=/3'=7=7'=l 



and so on. Note that the series representations of these parametric derivatives (obtained by 
differentiating the formula (|T3p) will contain ?/^-functions , 



az 



and their derivatives, because 



d 

— (a)j = {a)j {ipia + j) - i){a)) 



(16) 



(17) 



etc. The representation (p!^ ) gives us asymptotic expansion of J{s,t;m) (^) for large |x| and 
\y\ (with due regard for Inx and \ny terms). 

To be sure that the analytic continuation (0) corresponds to the same function (^, we may 
compare the results for x < 0, y < 0. To obtain "standard" representations for the functions 
on the r.h.s. of (111), it is convenient to use the well-known reduction formula, 



F2ia,P,P'; /3,7V,t^) = (l-u)~%Fi 



a, P' 

i 



l-u 



where 2-F1 is a usual Gauss hypergeometric function. Using the well-known parametric integral 
representation for 2F1, 



iF, 



a, P 
7 



r(7) 



TiP) r(7-/5) 



^ (1 - e^)" 



(19) 
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and its parametric derivatives, it is easy to obtain the following results: 

F2 (1/2, 1, 1; 1, 1| u, v) = ^ = , (20) 

daF2 = In/?,, , (21) 



A 



xy 



A ,„(^k±A) ,u,(|^)-L.(-|^)}. (23) 



Note that analogous results for d^'F2 and dad^'F2 can be obtained from (^) and ( [23| ) by 
substitution x ^ y. 

To calculate 9^(9^/ F2, it is impossible to use (ITBI), and we need to employ a double integral 



representation for F2 |T^. Differentiating it with respect to 7 and 7' and putting a = 1/2, /5 



/3' = 7 = 7' = 1, we get 



3 ' ' 



d^dyF2 = ^uv I I d^dri ln(l - ln(l - v) (1 - < - t;?])"^/^ . (24) 







Evaluating the integrals yields 



Inserting the expressions (pOD-(p3D and (^) into ([T^), we obtain the same result as in (^. 
This fact confirms the correctness of analytic continuation from the variables x and y to 1/x 
and l/y. 



4. In the present paper we have discussed different representations for Feynman loop dia- 
grams by taking, as an example, the box diagram contributing to the scattering of photons by 
photons. Hypergeometric representation (^) is convenient for studying asymptotic expansions 
and constructing analytic continuation to other variables (0). Standard representation (P) 
(involving dilogarithms of rather complicated arguments) is convenient in numerical calcula- 
tions. Anyway, the formula of connection (^) enables one to pass from one representation to 
another (when it is necessary). The formulae (^) and ( [T^ ) explicitly give all terms of asymptotic 
expansions for small and large values of |s| and respectively. Note that these formulae do 
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not cover all possible regions of Mandelstam variables (in other regions we need to use other 
analytic continuation formulae). For example, in the papers |jlO|, analytic continuation to 
physical thresholds has been studied. We consider the presented analytic continuation as an 
interesting example of using the multiple hypergeometric functions theory in loop diagrams 
calculation. 
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